Tan's universal contact and collective oscillations of strongly interacting Fermi gases 



Mark DelloStritto and Theja N. De Silva* 
Department of Physics, Applied Physics and Astronomy, 
The State University of New York at Binghamton, Binghamton, New York 13902, USA. 

We study strongly interacting two component Fermi gas near a Feshbach resonance. By using a 
ground state energy functional constructed based on asymptotic limits and Monte Carlo calculations, 
we calculate the contact, structure factor, and collective oscillation frequencies in the BCS-BEC 
^vq . crossover region. The calculated contact and structure factor show excellent agreements with recent 

experiments. We compare these results with a standard mean-field theory and find that the contact is 
proportional to the square of superfluid order parameter. Further, we present the chemical potential 
. and the polytropic index in terms of homogenous energy and the contact. 

I. INTRODUCTION 

As a result of easy controllability enabled by present day laser technology, experimentalists were able to achieve 
remarkable progress in trapping and cooling of atomic gases towards degenerate temperatures [lj]. The strongly 
interacting Fermi gas is one of the richest systems studied by the cold atom community. Although harmonically 
trapped quantum gases are very dilute systems, most of their properties are governed by the interaction between 
particles. At these characteristic densities at ultra-cold temperatures, only isotropic and short-range s-wave scattering 
between particles can take place. This scattering can be characterized by a single parameter, the s-wave scattering 
length a. Experimentally, the scattering length can be tuned by using magnetically tuned Feshbach resonance (FB). 
In the unitarity limit, where a is tuned to positive or negative infinity, the system is strongly interacting. At low 
temperatures and in the attractive interaction regime where a < 0, atoms form Bardeen, Cooper, and Schreifer (BCS) 
pairs so that the ground state is a BCS superfluid. In the repulsive interaction regime where a > 0, atomic potential 
supports two-body molecular bound state in vacuum. In this limit the ground state corresponds to a Bose-Einstein 
condensation (BEC) of these molecules. In between these two states, there is a smooth crossover when a changes its 
sign by passing through ±oo. At zero temperature, the only relevant dimensionless parameter for physical properties 
is £ = (fci?a) _1 , where kp is the Fermi wavevector. Because of the short range nature of the interaction between atoms 
in a dilute gas, the range of interatomic potential tq is much smaller than the average inter-particle spacing. 

The challenging problem concerning the system properties is how to find the ground state energy E(£). The problem 
C<") 1 does not have an analytical solution so one has to rely on numerical solutions such as Monte Carlo simulations. In 
general, the solutions at the limiting cases; such as a — > 0* can be obtained from the asymptotic behavior. At 
unitarity where £ = 0, a must drop out of the problem so there are only two energy scales available; Fermi energy 
€f = (3n 2 n) 2 / 3 H 2 /(2to) and the temperature T. Here m is the mass of the atom and n is the atom density in 
^vq ■ the system. Thus, dimensional analysis follows that the average single particle energy must have the form e = 
ei?/i(fcsT/eF)iWhere h(x) is a dimensionless function. The universal constant h(0) = 3(1 + /?)/5 has been calculated 
theoretically [H and is in good agreement with experiments Q. The universal many-body parameter ft — —0.56 for 
a unitary gas and /3 = for an ideal Fermi gas. The same argument holds for other thermodynamic quantities at 
t-H ■ unitarity. As the physical properties are independent of microscopic details, by studying a unitary Fermi gas, one 
learns about the equation of state of a generic strongly interacting gas, such as high density nuclear matter found in 
• ■ the center of neutron stars. 

There is, however, another form of "universality" valid not only at £ = 0, but also in the entire BCS-BEC crossover 
region. This universality originated from the fact that kpro -C 1, i.e, the range of the interaction is much smaller 
than the inter-particle distance. Using generalized functions to deal with the singularities associated with the zero- 
range potential, Shina Tan recently derived a number of universal relations for Fermi gases in the BCS-BEC crossover 
region Q. These exact universal relations (Tan relations) connect microscopic properties to thermodynamic quantities. 
The connection is made through a single quantity termed contact (C), which is defined as the high momentum (fc) 
tail of the momentum distribution n(k); 
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C = lim k 4 n(k). (1) 

k— >oo 

Remarkably, this contact carries with it all the properties and many-body physics of the system. The Tan relations 
are applicable for broad situations: homogeneous or trapped systems, many-body or few-body systems, superfluid 
or normal phases, and finite or zero temperatures. They were later rederived by using a renormalization scheme 
in the quantum field theoretical framework by using a lattice model to regularize the singularity iff], by using a 
nonlocal quantum field theory Q , and by using a Schrodinger formalism [1] . For a ground state of a 4-fermion system, 
the contact has been calculated numerically for different values of interactions and found that the Tan relations are 
consistent within the systematic errors [9( . The temperature dependence of the universal contact has been investigated 
by several authors using a high-temperature quantum virial expansion method [icj |. ab-initio results 11 [, and a 
combined mean field- numerical method (l^ |. Concurrent with the preparation of this manuscript ref. [l3| presented 
complementary BCS mean field results with the argument that the contact is the conjugate of the inverse scatte ring 
length. Very recently, the Tan relations were experimentally verified at University of Colorado, Boulder, USA 14 1 
and at Swinburne University of Technology, Melbourne, Australia [l5| • 

In this paper we investigate the strongly interacting Fermi gases at zero temperature focusing on the universal 
behavior. By using an energy functional e(£) constructed by Kim et al. [lfjj], first we calculate the contact in the 
entire BCS-BEC crossover region. Then we calculate the structure factor from the Tan relations. Further, we 
find that the contact calculated from usual mean field theory is proportional to the square of the superfluid order 
parameter and shows very good agreement with recent experimental data and the energy functional method. In 
addition, we generalize the theory to trapped gases using local density approximation and polytropic equation of 
state. In harmonically trapped systems, we calculate the contact, energies, and collective mode frequencies. We use 
both hydrodynamic theory and sum rule approach and provide detail discussion about the validity of sum rules for 
crossover Fermi systems. All the calculated physical quantities are presented in terms of the contact and homogenous 
energy. 

The paper is organized as follows. In the following section, we introduce the energy functional and calculate the 
contact for a homogenous gas as a function of £ . Then we present our calculation on structure factor. In section III, we 
use BCS mean field theory to calculate the contact. In section IV, we generalize our theory to a harmonically trapped 
Fermi gas. The section V is devoted to discuss the collective oscillation frequencies and derive various oscillation 
frequencies. Finally, in section V, we provide a discussion and our summary. 

II. THE ENERGY FUNCTIONAL AND THE CONTACT 

We consider a two component Fermi gas with both spin components equally populated at zero temperature. We 
neglect the interaction between same spin state, but the interaction between opposite spin states is considered in 
the form U(r) = {AnTi 2 a/m)8{r). The total number density n and the Fermi wavevector kp are related through 
kp = {i^n) 1 / 3, . The ground state energy per particle e(£) is determined by the dimensionless interaction parameter 
£. As £ changes from — oo to to +oo, the ground state smoothly changes from BCS superfluid state to BEC state. 

In the BCS limit a — > 0~, the energy per particle can be expanded in powers of l/£: 

e (« = 2 -(s + R + "") (2) 

In the unitarity limit a — > ±oo, the energy per particle can be expanded in powers of £ @, Q : 
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where the universal constant /3 is estimated to be /3 = —0.56 ;2|. In the BEC limit a — >• + , the energy per particle 
can be expanded in powers of £ -3 / 2 : 
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where a m ~ 0.6a is the boson-boson scattering length 17] . Here the boson is the diatomic molecule formed out of 
two fermions in the BEC regime. The leading term is the total binding energy density for molecules and the second 
term is the energy density of molecular BEC. Using the Monte Carlo results at unitarity and asymptotic results at 
BEC and BCS limits, the energy per particle in the entire BCS-BEC regime e(£) = £fM£ -1 ) nas been constructed 



from a Pade parametrization by Kim et al. [16|, [18j. In the BCS regime where £ < 0, 
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and in the BEC regime £ > 0, 



The coefficients are found to 
2.29188(0.200149), and <S 4 (a 4 ) = 



h{x) = _ 2 ^x+a^ (6) 
v ' 2e F 1 + a 3 x + a 4 x 2 y ' 

be 5i(ai) = 0.106103(0.0316621), 5 2 {a 2 ) = 0.187515(0.0111816), 5 3 (a 3 ) = 



1.11616(0.0423545). The molecular energy is E„ 



-h I mar. This functional 



has shown a very good agreement with constrained variational approximation [19| and correctly produces the limiting 
behavior. 

Adiabatic sweep theorem; one of the Tan's universal relations, connects the contact per unit volume (contact 
density) c = C/V and the energy of the system 
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The adiabatic sweep theorem gives connection between the change in total energy density due to adiabatic change 
in the scattering length and the contact density. As the contact measures the probability of two antiparallel spins 
being close together, the adiabatic theorem links the short-range behavior to thermodynamics [5j. By using the 
energy functional e(£), we calculate the contact density in the BCS-BEC region as a function of £. The calculated 



dimensionless contact density s = c/k F is shown in FIG[TJ Finite temperature experimental data from ref. 14 1 are 
shown in filled circles. While the black circles represent radio frequency measurements, the gray circles represent the 
momentum distribution measurements. Notice that the data in figure 2 of ref. 14j shows a different scaling because 



of their data represents the contact per particle c p — 3ir s. As can be seen, our zero temperature model gives an 
excellent agreement with finite temperature experimental data on the BCS side of the resonance. The experiment is 
done at very low temperature T = O.llTp, where Tp is the Fermi temperature. Our calculation indicates that the 
low temperature contact is not very sensitive to the temperature. The tiny bump at the BEC side of the unitarity 
in our calculation is an artifact originated from the Pade approximation used to construct the energy functional. 
Alternatively, one can use the energy functional proposed by Manini et al. [28[ to calculate the contact. Even though, 
Manini et al.'s energy functional is more accurate on the BEC side around £ = 0, it is discontinuous at £ = 0. 
Therefore, the contact shows a large jump at unitarity. In section III below, we compare the contact calculated in 
the present section with the BCS mean field theory. Later in section IV, we calculate the contact of harmonically 
trapped gas. 



Structure factor 

Tan also derived an universal relation for the large momentum behavior of the spin-antiparallel static structure 
factor S^(q) = S(q). This relation originated from the short-range behavior of the pair correlation function ni^ (r) = 

J dR(n-|-(R— r/2)n^(R+r/2)). Taking the limit r 0, and Fourier transforming of nffl(r) gives Tan's static structure 
factor 

S(^^(l-^4 (8) 
4 q \ tt q J 

This is the large q 3> &f behavior, however, q must be smaller than the inverse of range of the interaction potential. 
In FIG. [H we present the results for q = 5/cf in the BCS-BEC crossover region. The filled circles are the experimental 
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FIG. 1: The dimensionless contact density s = c/k% as a function of dimensionless interaction strength. The filled circles are 
experimental data from reference [lj]. Gray filled circles: measurements from momentum distribution. Black circles: radio 
frequency measurements. 




FIG. 2: The static structure factor S(q) as a function of dimensionless interaction strength. The filled circles are experimental 
data from references [l5|, H3] ■ The thin blue line is the theoretical calculation based on random phase approximation from 
references [HIH. 



data from Ref. [l5| and [20|. The thin line is the theoretical calculation based on a random phase approximation 
(RPA) from reference 21 1 and [2(J. Except at the BEC side of the unitarity, our model captures the universal power- 



law (1/q) tail of the static structure factor and shows an excellent agreement with experiments [l5[ and recent theory 
based on RPA method [211 ] . 



III. THE CONTACT FROM MEAN FIELD THEORY 



We restrict ourselves to the wide Feshbach resonance, so that we do not need to explicitly consider the closed 
channel molecule contribution. The fermions of different hyperfine states f an d i interact through a short-range 
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effective potential U(r) = gS(fl — r). The system of N atoms is then described by the Hamiltonian H = H\ + H2 
where H\ and H2 contain kinetic and interaction part of the Hamiltonian defined as, 



H x = E/ d 3 ri>Ur)l-^-»}Mr) (9) 

U 

H 2 = -g J d 3 rV4(f)^j(r)^(f)^(rO. 

The field operators ip<r(r) obey the usual fermionic anticommutation rules, and describe the annihilation of a fermion 
at position f in the hyperfine state a. The parameters m and [i are the mass and the chemical potential of the atoms. 



Using the usual BCS mean- field decoupling |22[, the zero temperature grand potential is given by 



(e k - M ) - V(tk - V? + A 2 
2 \2irfi 2 a V £ fc 



(10) 



The BCS-Bogoliubov excitation spectrum of the atoms is given by Ek = \J (ek — /-i) 2 + A 2 , where ek = h 2 k 2 /2m is 
the kinetic energy and A = g{ip^(r)i^f(r)) is the superfluid order parameter. Notice that we have already eliminated 
the ultraviolet divergences originated from the nature of short-range potential. This was done by regularizing the 
interaction term according to the relation 
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The gap equation is obtained by the minimization of the grand potential density with respect to the superfluid order 
parameter, 



The fermions number density n = N/V is obtained by the variation of grand potential density with respect to the 
chemical potential, 



ek~ A* 



y/fa - + A 2 

Finally, evaluating dA/da from the gap equation, we find the contact density 
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By converting the sum into an integral over the momentum, we numerically solve both gap equation and the 
number equation simultaneously to find the superfluid order parameter A. The mean field contact density is shown in 
FIGlH As can be seen, the mean field dimensionless contact shows a reasonable agreement with our energy functional 
calculation carried out in the previous section. The inset shows the comparison between BCS mean field theory 
and the experiments. At the BCS side of the unitarity [(kfa)^ 1 = —0.0386], the contact calculated from the energy 
functional method (0.0951) gives a remarkable agreement with the experimental value (0.0951). However, the BCS 
mean field result (0.1096) is off by 15% at (/c/a) _1 = —0.0386. Obviously, while the contact goes to zero in the weak 
coupling limit (a Q~), it approaches to infinity in the strong coupling limit (a — >• + ). 
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FIG. 3: The dimensionless contact density s = c/k% as a function of dimensionless interaction strength. The gray line is the 
BCS mean field results and black line is the results from energy functional method. The filled circles are experimental data 
from reference [3]. See FIG Q] also. 



IV. HARMONICALLY TRAPPED FERMIONS 

Using the short-range nature of the interaction, Tan derived an additional universal relation for fermions in a 
harmonic trap; 

h 2 

T -V + U = — C. (15) 

Here T, V, and U are the kinetic, harmonic trapping potential, and interaction energies respectively. This generalized 
virial theorem was first experimentally verified at unitarity by Thomas et al. 24] and later by Stewart et al |14| for a 
range of interaction strengths. In this section, we make use of generalized virial theorem in local density approximation 
to calculate the contact density in a trapped system. 



Chemical potential and polytropic index 

In terms of energy functional given above e(£) = eph, the chemical potential of the Fermi gas can be calculated by 
using Gibbs-Duhem relation 



(i = d[ne(£)]/dn. (16) 

Using this relation, we find the chemical potential /i = ep(5h/S + 7rs£/2), where s is the dimensionless contact density 
defined above. The calculated chemical potential as a function of interaction is shown in Fig|4j 

The polytropic index 7 relates the pressure and the density of the gas as P cx ri y+1 . Alternatively, 7 relates the 
chemical potential as fi cx n 7 . We calculate this effective polytropic index as the logarithmic derivative of the chemical 
potential 7 = (n/ (i)dfi/dn. In the BCS-BEC crossover region, we find that the effective polytropic index reads 



= l0h/3 + 3ir S £ + ir(d S /dZ- 1 )/2 
7 5h + 3irs£,/2 ' 1 ' 

As one expects, we find 7 ~ 2/3 at both non interacting limit and the unitarity limit 26]. As we have included 
the molecular energy in our formalism, we find 7 = in the deep BEC limit. Without inclusion of the molecular 
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contribution, we find 7 = 1 in the interaction dominated BEC regime. The effective polytropic index in BCS-BEC 
crossover region is plotted in FIG. [5] In the crossover region, the effective polytropic index which excludes the 
molecular binding energy varies between 7 ~ 0.6 — 1. 



Contact density in a trapped system 

In order to use the generalized virial theorem to find contact density in a harmonic trap, we first calculate the total 
trapped energy in the trap. Assuming the Fermi gas is locally homogenous, we use local density approximation and 
the effective polytropic index to calculate the trapped energy. Within local density approximation, the inhomogeneity 
is taken into account by using spatially varying chemical potential as [i = fi — muj^r 2 /2, where ui and /i are the 
trapping frequency and the chemical potential at the center of the trap. For example, the total number of atoms, 
N = J d 3 rn(r) can be converted into integration over \x using d[i — —muj^rdr. Assuming /1 = An 1 ', we find 




4\/27r Z 1 ^ 

N = w4^ Lj lhv ^ d »- (18) 

Similarly, the sum of the kinetic and interaction energies of the trapped system T+U = E = j d 3 re[£(r)] and trapping 
potential energy V = towq/2 J d 3 rr 2 n(r) can be converted into the integral over Evaluating the integral, we find 
E = (2/37) V. Combining all the energy components, we find the total trapped energy Et = T + U + V , 



The calculated total trapped energy in the BCS-BEC crossover region is shown in FIG. |U Finally, the ratio of the 
contact density to the interaction strength in tapped Fermi gas is calculated from the Tan's generalized virial theorem 
in Eq. JI5]), 



4£t 2-37 
37riYe F 2 + 37' 



(20) 



The contact density of a trapped Fermi gas is shown in FIG. [7] The filled circles are experimental data from 
reference 14j . 



V. COLLECTIVE OSCILLATIONS 



Hydrodynamic theory 



The dynamic of the gas can be described by the time-dependent nonlinear Schrodinger equation, known as the 
generalized Gross-Pitaevskii equation [25|, l29j 



V 2 + V(r)+/i[n(r,i)] 

2m 



(21) 
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FIG. 7: The dimensionless contact density in a trapped Fermi gas. The filled circles are experimental data from reference [lil ] 



where ip(r,t) = \J n(r, t) exp[i<p(r, t)] is the superfiuid wave function. In terms of superfiuid density n and velocity 
v = (7i/m)V^>, Eq. (|21j) can be converted into continuity equation and Euler equation, 



— = -V ■ [n(r)v] 



and 



<9v 

l 9F 
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V(r) + fj,(ri) + qp] 



(22) 



(23) 



where qp — (— ?i 2 V 2 v / n)/(2mi/n) is the quantum pressure term. In Thomas-Fermi approximation (TF), we assume 
that the gas is locally uniform so that the quantum pressure term can be neglected. This TF regime can be realized 
in a trap in the limit N — > oo. It is interesting to note that once we drop the quantum pressure term these equations 
are classical so that these hydrodynamic equations are applicable for both Fermi and Bose superfluids. However, the 
hydrodynamic equations are sensitive to quantum corrections, statistics, and dimensionality through equation of state 
which enters through the density dependent local chemical potential /x(n). This hydrodynamic description is valid as 
long as the collisional relaxation time is much smaller than the inverse of the collective oscillation frequencies. 

For the ground state n(r,t) — n (r) and v(r,t) = 0. In order to study the collective oscillations above the ground 
state, we linearized the hydrodynamic equations by substituting n(r, t) = n (r) + 5n(r,t) and v(r, t) — 5v(r,t). After 
neglecting the higher order terms and keeping only the linear terms and then by combining the continuity and Euler 
equations, one finds the linearized version of the hydrodynamic theory 



d 2 Sn 



,n dfi 
[ — V — dn\. 

m on 



(24) 



Taking 5n(r, t) = Sn(r) exp[zu;i] and u = An 1 ', this equation can be converted into a linear eigenvalue problem. In the 
TF regime, the continuity and Euler equations admit the ground state density no(r) = no(l — r 2 / R^p) 1 / 1 , where the 
Thomas- Fermi radius Rtf = (mwg). Re-scaling the length, density, and oscillation frequency by introducing 

f = r/RxF, no(f) = no(r)/no(0), and uj = lu/ujq respectively, Eq. (|24|) can be written in dimensionless form, 



2u 2 



Sn = V ■ [noVnQ Sn] 



(25) 
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FIG. 8: The monopole mode frequency of a Fermi gas in a spherically symmetric trap. At weakly interacting BCS limit 
lom — 2cj and deep BEC limit lom = v5w . At unitarity, mode frequencies are equal to the ideal gas limit required by the 
universality. 



For a spherically symmetric trap, the angular momentum I and its projection in the z-axis m are good quantum num- 
bers. Following the Ref. (2|| and Ref. (27|, the solution for the density fluctuation have the form Sn — R(r)Yi m (0, </>), 
where = f i nj! ) _7 (l — f) 1 / 7-1 f(f). The function /(f) satisfies the hypergeometric differential equation, 



f(l 



f) — - + 

' dr 2 



Z + 3/2 - (Z + 3/2 + l/ 7 )f 



+ - 



w 2 



'2; 



dr 

l -f. 



(26) 



The solutions of this standard equation is hypergeometric functions [26l . 1271 . l29j and the eigenvalues are given by 



uil l = l+n r j(2n r + 2l + l) + 2n r , (27) 

with radial quantum number n r — 0,1,2,... and angular momentum quantum number I = 0,1,2,... Using the 
polytropic index calculated in the previous section, all the modes of collective oscillation frequencies can be presented 
in terms of the contact density from this energy eigenvalues. For example, monopole mode frequency or the breathing 
mode frequency is given when n r = 1 and I = 0. Zero temperature monopole mode frequency ujm = wo\/37 + 2 of the 
two component Fermi gas in a spherically symmetric trap within the BCS-BEC crossover region is shown in FIG. [8] 
The frequencies are given as a function of the dimensionless interaction parameter £. Our model predicts ujm = 2lo d 
at the weak coupling BCS limit. At the weakly repulsive BEC limit, ujm — y5oj . These results are consistent with 
the sum rule approach present in the next subsection. At unitarity, ljm = 2u as required by universality. The other 
mode frequencies, such as second breathing mode (n r = 2, / = 0) and quadrupole modes (/ = 2) can be readily 
calculated from Eq. (|27l) . 



Sum rule approach 

In this section, we use the sum rules approach in the linear response theory to express rigorous upper bounds to 
the energies of the collective oscillations at zero temperature. The response function xi z ) a t complex frequency z is 
written in terms of the imaginary part of the response function x"( z ) as 
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f duj y"(z) 

X (z) = / —^4- 28 

J 7T U) — Z 

Expanding the response function in terms of 1/z, the moment expansion of x( z ) is given by 

= E \ m v-^ ( 29 ) 



where the pth moment is defined as 



zP 
p 



— ^x»- (30) 

7T 



In general, one has to calculate the response function \ associated with the mode excitation operator to evaluate the 
moments. However, these moments can be computed without calculating the response function if one uses the sum 
rule approach [23J. A major advantage of sum rule approach is that the collective mode frequencies can be calculated 
without the full solution of the eigenstates of the Hamiltonian. One disadvantage is that this method provides the 
rigorous upper bounds for the lowest collective mode frequency excited by the excitation operator. The excitation 
energies of the collective modes are given by the ratio m p ^-i/m p or y / TO p+ 2/m p . When the response function is almost 
exhausted by a single mode (that is, the state is highly collective), then these ratio giv e the exact collective mode 
frequency. We estimate the excitation energies of the collective modes from the ratio [32j 




Using the completeness relations, the moments can be written as 

mi - \{[Q*,[B,Q]]) 

m 3 = 1([[Q\H],[[H,[H,Q]]]) (32) 

where [A, B] represents the commutator between the operators A and B, Q is the mode excitation operator, and H 
is the Hamiltonian of the system 



= ^(T a + V a ) + U (33) 

ot 

where a = x, y, z. The excitation operators corresponding to the dipole, monopole and quadrupole modes are 
Q = J2i z i> Q — J2i r i> an d Q — J2i( r i ~ 3zf) respectively. By evaluating the commutators above, the collective 
mode frequencies for the dipole, monopole, and quadrupole modes are [13] 



wd = u (34) 

_ It + w + w 

- 9 / T + V 



Notice that the dipole mode is not effected by the interaction. This mode corresponds to the oscillation of the center 
of mass of the atomic cloud. The calculation of monopole collective mode frequency from the first and third sum rules 
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involves the evaluation of interaction energy U of the trapped system. The calculation of this interaction energy for 
a crossover Fermi gas, especially at unitarity is not trivial. Indeed, the use of third sum rule for crossover fermions 
is problematic. For example, consider the equation of motion of the operator Q = Q a = X)i r L- Using Heisenberg 
equation of motion, one finds Q a = [Q a , H]/(ih) = Y^i( r iaPia +Pia r ia)/ m an d Q a = 4(T Q — V a + U/2)/m. For 
a spherically harmonic oscillator potential, the three terms represent the average kinetic energy T — (X^Pi /(^ m ))i 
the harmonic oscillator potential energy V = muj 2 J r 2 n(r)d 3 r /2, and the average mean-field interaction energy 
U = g J n(r) 2 d 3 r /A. In equilibrium Q a is time independent, so Q a — 0. Then summing over all three components, 
this follows the original virial theorem, 



2T-2V + 3U = 0. (35) 

Now, let's combine this original virial theorem in Eq. Q35[) with the Tan's generalized virial theorem given in Eq. (|15l) . 
For both non-interacting fermions and unitary fermions, the interaction energy U turns out to be zero at unitarity. 
In other words, within our sum rule approach, the third sum rule (together with original virial theorem) leads to an 
incorrect interaction energy at unitarity. Regardless to this incorrect interaction energy, the combination of these two 
virial theorems [Eqs. (|15[) and (1351) ] into the sum rules produces the correct breathing mode frequency. In any case, 
we avoid using Eq. (1351) , instead we use compressibility sum rule to evaluate the breathing mode frequency. Contrast 
to the other sum rules, the compressibility sum rule 



m_ x = lim f J^Le^x"^), (36) 

cannot be written as a set of commutators. By completing the Fourier transform over the frequency (or using the 
Kramers-Kronig relation), one finds that the compressibility sum rule m_i = x(w = 0) = dQ/db is the static response 
function, where b — mwg/2. Using ljm = \J m\/m,-\ for the breathing mode frequency, we find 



_ -A(r 2 ) 
ma(r^} /db 



Defining (r 2 ) = An r 4 n(r)dr, where ro = y/no/b and using /i = An 1 , we find 



d(r 2 ) = 2(r 2 ) 
db 6(3 7 + 2)' 1 ' 

The lowest order breathing mode frequency for spherically symmetric trap is then given by % = + 2. This 

is exactly the same results predicted by the hydrodynamic theory in previous subsection. This indicates that the 
rigorous upper bounds for crossover Fermi gases calculated from the sum rules are very tight. In other words, the 
lowest order breathing mode frequency of a trapped Fermi gas is highly collective. 



Collective modes of a Fermi gas in an axially symmetric trap 



In an axially symmetric trap with trapping frequency V(r) — m[uj r (x 2 + y 2 ) + u> z z 2 )]/2, the monopole modes 
frequencies can be calculated easily using the polytropic index 7. Generalizing the hydrodynamic theory for an axially 
symmetric trap, the lowest frequency axial monopole mode and the radial monopole mode cu a = lj z 
T) respectively [H,^ 

n 2£ 



The other limit where A > 1, the monopole modes are ui a 



y/3- (7 + I)" 1 and 

These expressions are valid only for a cloud with aspect ratio A = uj z /uj r <C 1. 

3 



uj z J2T^ and cj+ = w r% /(67 + 4) (7 + 2) 
The experimental results are available only for the case A< 1 (IS |l| . These two mode frequencies in the BCS-BEC 
crossover region are shown in FIG. [9] For the limiting cases when 7 = 1 and 7 = 2/3, the axial mode frequencies are 
uj a = co z ^5/2 and u> a = lo z ^12/5 and the radial mode frequencies are w+ = 2cu r and uj a — cu r a/10/3 [32j |. 
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FIG. 9: The axial (a) and radial (b) monopole mode frequencies in the BCS-BEC crossover region. 



VI. SUMMARY 



By using two methods; a ground state energy functional constructed based on asymptotic limits and Monte Carlo 
calculations and the BCS mean field theory, we calculate the contact of a two component Fermi gas near a Feshbach 
resonance. The calculated contact which encapsulates all the many body physics in the BCS-BEC crossover region 
shows excellent agreement with the recent experiments. Within the mean field theory we find that the contact is 
proportional to the square of the superfluid order parameter. We then investigate the large momentum behavior of 
the structure factor using a Tan's universal relation and the contact. Our structure factor calculation in the BCS-BEC 
crossover region also shows an excellent agreement with recent experiments and a RPA based theoretical calculation. 

We combine the Tan's generalized virial theorem and local density approximation to derive the contact density in 
harmonically trapped fermions. Using the polytropic equation of state, we calculated the trapped energies and the 
collective mode frequencies of a Fermi gas trapped in a harmonic oscillator potential. We present all these quantities 
in terms of the homogenous energy density and the homogenous contact density. 
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